The Physics of Chromospheric Plasmas 
ASP Conference Series, Vol. 368, 2007 
Petr Hemzel, Ivan Dorotovic and Robert J. Rutten, eds. 



The Fine Structure of Solar Prominences 



Petr Heinzel 

Astronomical Institute AS, Ondfejov, Czech Republic 

Abstract. Solar prominences and filaments (prominences projected against 
the solar disk) exhibit a large variety of fine structures which are well observed 
down to the resolution limit of ground-based telescopes. We describe the mor- 
phological aspects of these fine structures which basically depend on the type 
of a prominence (quiescent or active- region) . Then we review current theoret- 
ical scenarios which are aimed at explaining the nature of these structures. In 
particular we discuss in detail the relative roles of magnetic pressure and gas 
pressure (i.e., the value of the plasma-/?), as well as the dynamical aspects of the 
fine structures. Special attention is paid to recent numerical simulations which 
include a complex magnetic topology, energy balance (heating and cooling pro- 
cesses), as well as the multidimensional radiative transfer. Finally, we also show 
how new ground-based and space observations can reveal various physical as- 
pects of the fine structures including their prominence-corona transition regions 
in relation to the orientation of the magnetic field. 



1. Introduction 

Solar prominences have been observed since the invention of the spectrohelio- 
scope and as an example one can see we ll-known Ha drawings by iSecchi. (18771 ) 



(one can be found in the textbook of Tandberg-HanssenI ( 19951 )). Although 



these pioneering observations were rather simple, they already indicate that 
prominences consist of many complicated structures having seemingly chaotic 
behaviour. Prominence fine structures have then been frequently observed with 
better and better resolution, reaching today a hundred of km. However, although 
the fine structure was known, the prominences have been modeled for decades 
as a whole, using simpHfied one-dimensional (I D) slab models. Their magnetq - 
hydrostatic (MHS) structure was first derived bv lKippenhahn Sz Schliiter ( 19571 ) 



(hereafter referred to as KS-model) and their radiative properties were stud- 
ied in a ID-slab approximation by many authors. Such models reproduced 
rather well the spectral properties of prominences as observed with lower res- 
olution. The aim of this review is to discuss various aspects of current inves- 
tig ations of prom ii ience fine st ructures. This topic was thoroughly reviewed 
bv lHeinzel Vial ( 19921 ). while Engvold ( 20041 ) summarized the latest observa- 



tional results achieved with the highest resolution. Some aspects of prominence 
fine structures related to space research were briefly discussed by IVial (2006). 
For a more general d escription of the promine nce physics the reader is referred 
to the monograph of Tandberg-Hanssen ( 19951 ). 
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2. Morphology of Prominence Fine Structures 

Prominence fine-structure morphology manifests itself rather differently in case 
of the limb observations and in case of disk filaments. Moreover, one has to 
clearly distinguish between typical quiescent and active-region prominences or 
filaments. Generally speaking, a prominence seen on the limb has appeared 
before or will appear later as a filament on the disk. In case of quiescent promi- 
nences larger-scale structures remain fixed while the fine structure changes rather 
rapidly. However, it seems to be very difficult to identify the same structures 
as seen on the limb and on the disk - this makes a lot of confusion and we will 
discuss this later. 

Our information about the fine structure morphology and dynamics comes 
from currently available high-resolution ground-based observations mostly made 
in the Ha line. With the instruments like the new SST (Swedish 1-m Solar Tele- 
scope) or DOT (Dutch Open Telescope), one can see fine structures down to 
the resolution limit (0.15" for SST or around 100 km). Homogeneous time se- 
ries are now expected with a similar spatial resolution from the Solar Optical 
Telescope (SOT) onboard the Hinode satellite. A large variety of fine structures 
and their dynamics is also seen on TRACE movies, although the spatial resolu- 
tion is lower, around 1" (see e.g., TRACE filament movies on DVD provided by 
LMSAL). These images are usually taken with a 171 or 195 A filter where the 
hot coronal structures appear simultaneously with cool ones - see Fig. [TJ Cool 
prominences or filaments are dark against the bright background which is due 
to the absorption of the background coronal radiation emitte d in these lines by 



the hydrogen and helium resonance continua (see cartoon in iRutt en 1999) and 
partially due to lack of emissivity of the TRACE lines within a volume occupied 
by cool prominence plasmas. At 171 or 195 A, the He I and Hell absorption 
dominates over the HI and it was s hown theoretically tha t this opacity is quite 
comparable to that of the Ha line ( Anzer Sz Heinzel 20051 ). 



A large quiescent prominence observed at Big Bear Solar Observatory (see 
Fig. [2]) exhibits many vertically orie nted threads of the cool plasma which are a 
few hundred kilometers wide (from Low &: Petrid 20051 ). On a larger scale, one 



can identify a few vertical plasma sheets which, using a lower spatial resolution, 
would appear as more-or-less homogeneous slabs with the thickness typically 
smaller compared to the other two spatial dimensions. Such a kind of low- 
resolution images led modelers to use a vertical ID slab approximation to a 
real prominence geometry which completely neglected the fine structure. In 
case of MHS models we already mentioned the classical KS solution. To handle 
the non-LTE radiative trans f er, ID slab models were ext ensively used starting 
from [Poland fc Anzeil (|l97ll ). lYakovkin fc Zel'dinal (|l975l ) and others. We will 



not detail these models here but rather discuss their recent generalizations to 
fine-structure modeling. 

On the disk, the high-resolution Ha images or movies show fine-structure 
fibrils of different lengths, the thinnest visible down to the resolution limit of SST 
or DOT. Very thin dark fibrils (we call them "fibrils" to distinguish them from 
vertical "threads" seen on the limb) visible along the spine of a quiescent filament 
are rather short and inclined to the filament axis due to the shear of the magnetic 
field lines (Fig. [3]). On the other hand, much longer fibrils can be seen within the 
barbs or connecting various parts of the filament main body (Fig. U]). The fact 
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Figure 1. TRACE image of a limb prominence taken in the 195 A line. 
Note the fine structures which appear very dynamic on the respective TRACE 
movies. 

that densely packed fibrils seen along the spine are rather short compared to a 
large-scale magnetic arcade in which the prominence/filament sits indicates that 
these fibrils are locations of cool plasma condensations in a presumably dipped 
magnetic field. 



3. Dynamics of Prominence Fine Structures 



While the large-scale quiescent prominence structure is rather stable, the fine 
structures described above exhibit a strongly dynamical behaviour. Their shape 
and brightness change on time scales of minutes which was already noticed by 
Engvoldl (11976') using the Ha prominence observations made at the Dunn VTT 



at Sacramento Peak. Today one can study the fine-structure dynamics and 
prominence evolution on prominence/filament movies taken by TRACE in the 
195 A line or on high-resolution Ha images or movies. Within the disk filaments, 
individual fibrils move sideways with vel ocities up to 3 km s~^ which se ems 



to be consistent with limb observations of lZirker Sz Koutchmvl (i99(]| . 199ll ). A 
qualitatively new observation was reported bv lZirker et al.l (|l998l ) who identified 
a kind of streaming and counter-streaming in a filament in both spines and barbs. 
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Figure 2. Big Bear Solar Observatory (BBSO) image of a typical qui- 
escent prominence . Quasi-vertical plasma threads are well visible (from 
ILow fc Petri"i[200l . 




Figure 3. Dutch Open Telescope (DOT) image of a disk filament. Many 
dark fibrils are seen along the filament spine and elsewhere (courtesy of R.J. 
Rutten). 



having flow velocities around 15 it 10 kms ^. 
for distance s of 10 ^ to 10^ km at these speeds 
Zirker et al.l (|l998l ). 



Dark Ha knots were tracked 
- see long arrows in Fig. 2 of 
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Figure 4. Swedish 1-m Solar Telescope (SST) image of a disk filament. Here 
we can see long dark fibrils which resemble rather thin flux-tubes (courtesy 
of O. Engvold). 



Ha 3VST. [9m-00-L9 07;aV;4GUT rrace:7lA 1998-06-19 a7'38:aiUT 




Figure 5. Ha and TRACE fine-structure fibrils. The spatial resolution of 
TRA CE is lower, but s ame dark fibrils can be well identified on both images 
(from ISchmieder eFall ^004) . 



The same counter-streaming was further observed bv lLin et ah ( 2003) us- 



i ng th e hig h-resolutio r i SST Ha images and Dopplergrams. Both IZirker et al 



( 19981 ) and iLin et al.l ( 20031 ) interpret these motions as plasma flows along the 
magnetic flux tubes which are not dipped. We have an evidence that these 
high-resolution Ha images show fine structures which are well coaligned with 
dark features visible in the s ame filament on TRACE 195 A images - see Fig. [5] 
from Schmieder et al. ( 20041 ) . These authors have studied the same filament as 



Lin et al.l (|2003l ). 
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4. Thermodynamic Properties from Spectral Diagnostics 

Spectral diagnostics of prominence fine structures is a difficult task for several 
reasons. First, one has to distinguish between optically-thin and optically-thick 
fine structures. In the former case several fine-structure elements (FSE) are seen 
along the line of sight and their radiation output is thus integrated. On the other 
hand, in case of thick structures, we mostly see only one FSE, most probably 
located closer to the border of a prominence rather than in its central parts. The 
situation is actually more complex, because an optically-thick FSE is thick in the 
line core but becomes thin in the line wings, where we again see many FSE along 
the line of sight. This situation largely complicates the analysis of spectral data. 
Second reason is that, as already mentioned above, FSE are highly dynamical 
and thus any multi-wavelength observation has to be made simultaneously or 
at least quasi-simultaneously (this means within time shorter than the life-time 
of the FSE) in all wavelength bands of interest. Another complication arises 
from the lack of our detailed knowledge about the internal structure of FSE 
and their distributions within a prominence - even their elementary geometry 
is subject of a controversial debate (vertical threads, horizotal fibrils with flows, 
blobs etc.). For spectroscopic work, two basic models are usually considered: 
(i) Each FSE has its own prominence-corona transition region (PCTR), its 3D 
shape will then depend on the magnetic field threading the FSE because of prin- 
cipally different thermal conductivity along and across the field lines, (ii) FSE 
are more-or-less isothermal, but their temperature increases towards the promi- 
nence boundaries and this forms another kind of PCTR on a larger scale. A 
combination o f (i) and (ii) is also possible. For further details and references see 
the review bv IVial ( 19981 ) and his recent summary at the SOHO-17 workshop 



(|via]l2006l ^. The scenario (i) was co rroborated by various author s , whi le the 
second one (ii) was suggested e.g., bv iPoland Sz Tandberg-HanssenI ([1983) who 
analyzed the spectral data from SMM/UVSP. Here it is worth mentioning that 
the scenario (ii) allows one to use a I D slab model with a re levant PCTR on 
each side as a first approximation (see lAnzer fc Heinze]|[l999l '). while the effect 
of a fine-structure PCTR (case (i)) has to be modeled by considering a super- 
position of several FSE (each of them, however, can be again mo deled by a 
simple ID slab w'hich desc ribes, contrary t o (ii) , only one FSE) ~ see lVial et al.l 
(|1989I 1. iHeinzell (|l989l ^ or iFontenla et al.l (|l996l ). Finally a principal difficulty 



in our understanding of the thermodynamical properties of prominence plasmas 
comes from the nature of radiation processes involved. The prominence spectral 
lines are mostly formed, namely in cooler parts, by the scattering of the inci- 
dent radiation coming from the solar surface (photosphere, chromosphere, or 
even PCTR, depending on the spectral line under consideration). This leads to 
strong departures from LTE and thus the non-LTE theory of radiative transfer 
must be applied to prominences and filan ients. This is known fo r a long time, 
here we can mention pioneering works b v iPoland Anzed (Il97ll 1. by an HAO 



group in the s eventies (see.H easlev & Mil kevi (119781') arid refe rences therein) or by 



Kiev groups (lYakovkin fc Zerdinalll975l . lMorozhenkolll984l) . For a recent sum- 
mary of modern non-LTE techniques see iHeinzel &: Anzed (|2005l l. The use of 



these transfer methods for ID slab s is almost routine today and s everal ID non- 
LTE codes do exist : the IAS code ( Gouttebroze Sz Labrossell200Cll ) . the Ondfejov 



code (lHeinzellll995l ). PANDORA is now also being modified to prominence slabs 
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1982 ) 



19941 ) 



(|Heinzel AvrettI 120071 ') 2D slab mod els were conside red first by IViaM 

usin g the M AM code of lMihalas et al.1 (IT978I). then by lAuer k Paletoul . 

and IPaletoul (|l995l ) and finally by lHeinzel Anzerl (|200lh . All these codes can 
be or haye been used to model either the prominence or filament as a whole 
(note: for prominences one uses slabs standing yertically aboye the solar sur- 
face, while in case of filaments the slabs are horizontal, i.e., parallel to the solar 
surface) or indiyidual FSE. Howeyer, in case of seyeral FSE, the situation be- 
comes much more complex beca use of thei r mutual radi a,tiye inte raction. This 
was treated in simplified ways by lZharkoyal ([1989 ) and byH einzell (|1989, ). but in 
most other cases the radiatiye interaction was neglected, all FSE were considered 
the be identical (as ID slabs) and simply superimposed along the line of sight 
( Fontenla et al.|[l996l : iHeinzel et alll200ll ). This interaction could, howeyer, be 
tackled using todays high-performance parallel machines. 

In the aboye summary we haye tried to giye the reader an idea how com- 
plex the spectral diagnostics of prominence fine structure is. Howeyer, we still 
completely neglected the dynamics or temporal yariations of thermodynamic 
properties. The actual yalues of thermodynamic parameters like the kinetic 
temperature (T), gas pressure (p), plasma density (p), electron density (ng), 
but also the number of FSE along the line of sight or their geometrical dimen- 
sions, deriyed froni promi n ence / f ilame nt spectra, haye been summarized e.g., by 
iTandberg-Hansse nI (Il995l). IViall (|l998l ) or recently, using SOHO/SUMER data, 
by Patsourakos fc Viall (|20oi r One approach is to study in detail the proper- 
ties at a giyen prominence location, where we see - as discussed aboye - one 
or more FSE along the line of sight, or to perform a statistical analysis of the 
raster data to get the information about a globa l distribution of su ch param- 
eters. Here we can meiition t he recent work by IWiehr et al.l ( 20071 ) (see also 
iLopez Ariste &i Aulanieil 12007 on p. 291 ff in this yolume), where the line ratio 
technique (HeDa oyer H/3) was used. The results based on high-resolution filtro- 
grams taken on VTT at Tenerife do indicate distinct situations: The case when 
for example the gas pres sure would be alm ost constant oyer the prominence 
region (see also iStellmac her fc Wiehr ( 20051 )) is of a particular interest since it 
indicates an almost linear decrease of the magnetic field intensity with height 
and a statistical homogeneity of the spatial distribution of FSE. 



4.1. Spectral diagnostics with SOHO 



A few months before this Coimbra meeting solar physicists celebrated 10 years 
of successful obseryations with the Solar and Heliospheric Obseryatory (SOHO), 
which is a joint ESA/NASA mission (see proceedings from SOHO 17 Work- 
shop held in Sicily, May 2006 - ESA SP-617). During this period, SOHO was 
frequently used to obserye solar prominences and filaments in UV and EUV 
spectral bands. In particular SUMER and CDS spectrometers proyided us with 
a wealth o f unique data. The use of SUMER for prominence physics was re- 
yiewed by iPatsourakos Vial ( 2002 ) . prominence spectroscopy with SUMER 

Con cernin g the prominence fine 



was also s ummarized bylHeinzel et al 



promm 
(12001) ■ 



structure, IWiik et all (|l999l ) and Cirigliano et al. (2004') haye studied the be 



hayiour of lines formed in PCTR, while iHeinzel et al. (2006) giye a summary 
of their work on the hy drogen Lyman spectrum. Th e Lyman lines of hydro- 
gen were also studied by lAznar Cuadrado et al.l (|2003l ). Finally, let us mention 
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that iParenti et al.l ( 20041 ) and iParenti et alj (|2005l ) have produced an atlas of 
the prominence EUV spectrum observed by SOHO/SUMER. 



5. Magnetic Field Determinations 

Magnetic fields in solar prominences were measured since several decades, using 
first the Zeeman technique and later on th e Hanle effect. This is summarized 
in the review by iLopez Ariste &: Aulaniei ( 20071 ) on p. 291 ff in this volume. 
Older polarimetric observations made with the coronagraph led, in the case of 
quiescent prominences, to magnetic fields typically below 20 G. For latitudes 
larger than 35° one obtains a mean value of B equal to 8 G , while at lower 



latitudes the mean field is around 10 G ( Bommier et al. 19941 ). Lerov et al 



( 19831 ) analyzed data for a l arge number o f pola r-crown filaments and found B 
in the range 2 - 15 G, while lAthav et al.l (|l983l ) arrived at values 6 - 27 G for 
a sample of 13 prominences (however, they didn't distinguish between quiescent 
and eruptive prominences). Practically all these measurements have indicated 
the predominance of horizontal fields which naturally posses the question why 
the fine-structure threads appear quasi- vertica l and not aligned along the field 
lines like the magnetic flux-tubes ( Lerov 19891 ). 

More recent meas urements seem to in dicate higher fields, reaching values 
of a few tens of Gauss (jLopez Ariste &: Aula nier 2007). However, the respective 
observations concern only a rather restricted sample of prominences which may 
not always be of a quiescent type or are restricted to low parts of prominences 
(including their feets). Higher fields also follow from the analysis of the Stokes 
V - signal which was not used in earlier studies. The actual intensity of the 
magnetic field is critical for our understanding of the prominence support and 
we will discuss this in the next Section. 



6. Magnetic Models of Fine Structures 

Since the fine-structure magnetic fields have not yet been detected (not because 
of lack of the spatial resolution, but due to still a low signal-to-noise ratio in 
polarimetric signals), most of our current knowledge or ideas come from data- 
driven modeling. Either one can use the detailed distribution of surface magnetic 
fields (like the SOHO/MDI maps) and p erform coronal extrapolat i ons of various 
degrees of sophistication (see review bv iLopez Ariste &: Aulanier ( 20071 ) in this 
volume), or use indirect spectroscopic methods and modeling to in f er the fine 
structure of the magnetic field (l i ke tha t suggested in iHeinzel et al.l (j200ll ) and 



recently used in Schmieder et al. ( 20071 )). Here we are mainly interested in the 



topolog y of the field i nside and in the surroundings of FSE (for a recent summary 



see 



also|A 



nzer 



(|200l )). 

Two different magnetic configurations are consid ered in relation to promi- 
nence fine structures. The first one, also discussed bv lLopez Ariste &: Aulanier 
( 20071 ). is based on the assumption that the plasma-/? (ratio of gas to magnetic 
pressure) is always very low, of the order of 10~^ or lower. In that case, the 
plasma itself has negligible effect on the topology of magnetic fields related to 
FSE. If, however, the plasma-/? is larger, of the order of 0.1 - 1.0 or even higher. 
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the weight of the plasma will produce what we call "gravity-ind uced magnet ic 
dips". Let us discuss these two situations in more detail (see also lAnzeill2002l ) . 



6.1. Low-/? models 

The to pology of the prominence magnetic fields at low /? is discussed by lLopez Ariste Sz Aulaniei 



( 20071 ) on p. 291 ff so that we give only a brief overview here. Today var- 



ious authors try to model the prominence magnetic field by linear or non- 
linear force- free-field ( "fff " ) extrapolations of the measured photospheric field. 
Aulanier &: DemoulinI ( 19981 ) show that the prominence may consist of many 



rather shallow magnetic dips which can be filled by the cool material visible 
e.g., in the Ha line. The p resence of such dips was further dem onstrated by 
Aulanier &: Schmieder ( 20021 ) and by I Aulanier &: Demoulin ( 20031 ). In the latter 



paper the authors studied several prominences and found that the extrapolated 
field strength is rather low for typically quiescent prominence (less that 10 G), 
higher than 10 G for a plage filament and can reach 30-40 G for active-region 
prominences. This clearly shows how important the knowledge of the promi- 
nence type is: when we speak about the fine structure of quiescent prominences 
(as in this review), we cannot use measurements obtained for other types of 
structures. Some of the recent measurements which indicate rather strong fields 
seem to be related to more active structures. Moreover, these new measurements 
were not obtained using the coronagraphs and thus, because of scattered light 
higher above the limb, the measurements are more restricted to lower heights 
where the prominences appear brighter and where the polarization signal/noise 
ratio is sufficient (Lopez Ariste - private communication). The usefulness of 
the filament m odeling with linear fff was rather convincingly demonstrated by 



Aulanier et al.. (2000) who performed a kind of "blind test": magnetic dips ob- 
tained from numerical extrapolation were marked by black bars schematically 
indicating the location of the absorbing Ha material ~ see Fig. O The shape of 
the filament and its various parts were then compared with the true shape as 
observed on the disk in Ha and t he agreement is quite re asonable. Non-linear 
fff models were considered e.g., bv lvan Ballegooiie n' (2004')^ 



For a given curvature of the magnetic dip, Aulanier Sz Schmieder ( 20021 ) 



assumed that the dip is filled by cool plasma up to one pressure scale-height H 
kT/gfi, where fi is the mean molecular mass (Fig. [7]). For typical prominence 
conditions H is of the order of 200 km. With the computed curvature of the fff 
dips one can draw the bars which indicate the presence of the absorbing material 
and have the length 21 as indicated in Fig. [71 However, the real distribution 
of the opacity will depend on the actual dip configuration, mass loading in it 
and on various plasma parameters including the hydrogen excitation/ionization 
conditions. This was not so far considered in the context of fff models. 

6.2. Gravity-induced dip models 

Spectroscopic determinations of the plasma density and its ionization degree, 
together with spectro-polarimetric determinations of the magnetic field lead in 
many cases to rather high values of the plas ma- /3. This is fo r exam ple the case 
of a sample of 14 prominences studied by iBommier et al.l ( 19941 ). where the 



plasma-/? reach es unity in severa l cases (see also the non-LTE analysis of these 
prominences bv lHeinzel fc Anzerl (|l998l ')^. Using a simple ID MHS configuration 
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Figure 6. Computed dark-fibril pattern from linear fff extrapolations (left) 
and the corresponding filament observed i n the Hct line (right). See the text 
where this "blind" test is described. From lAulanier et al.l ()200Cl[ ). 
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Figure 7. A schematic model showing how the magnetic dip is filled by 
the plasma up to one pressure scale-heig ht H. A cor responding horizontal 
extension of the mass loading is 21. From lPries^ (|1990f ). 



for a vertical fine-structure slab (representing vertical threads as seen in Fig. [2]) , 
one can easily derive a relation between the plasma-/? and the angle V' of th e 
magnetic-field inclination at the border of the dip (e.g., iHeinzel &: Anzer 19991 ) . 
P ~ cot^ 4^. 

Construction of equilibrium configurations in which the weight of the plasma 
inside the magnetic dip is balanced against the solar gravity by the Lorentz 
force is a difficult task. An analytical solution to this problem, assuming that 
the whole prominence is represented by a ID vertical plasma slab hanging in 
a dipped magnetic field, was proposed already in 1957 by Kippenha h n and 
Schliiter (so-called "normal" polarity model), while iKuperus &: Raadul ( 19741 ) 
have su ggested a similar model b ut with an "inverse" polarity. Based on the 
idea of iPoland &: Mariskal ( 19881 ) that the fine-s tructure threads a.re iii fact 
vertically- aligned magnetic dips loaded with mass, lAnzer fc Heinzell (|l999l ) and 
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Heinzel &: Anzej (|200ll ) started to model such gravity-induced dips using gen- 
eral ID or 2D MHS equilibria, respectively. Such local solutions within the 
magnetic dip are sometimes called KS-type models (according to Kippenhahn 
and Schliiter) because they use the KS-type analytical solution of the pressure 
equilibrium (see Section 8). However, they have nothing to do with the global 
magnetic topology of the whole prominence so that the inverse-polarity promi- 
nences can be treated in the same way. A comprehensive d escription of the 
ravity-induced dip models is given in the lecture notes by iHeinzel &: Anzer 
2005 ) to which the reader is referred for further details. Here we will only men- 
tion one important result, namely the dependence of the plasma-/? on the mass 
loading M and the fi eld strength B, which can be derived for ID vertical slab 
in MHS equilibrium (| Anzer fc Heinzell [20071 ) 



(1) 



Using the results of iGouttebroze et al.l ( 19931 ) . one can relate this M to optical 
parameters li ke the Ha line- c enter optical thickness or Ha integrated intensity. 

Recently, ILow Sz Petrid ( 20051 ) have also considered a series of ID gravity- 
induced dips aligned along the prominence spine and being in MHS equilibrium 
of the KS-type (Fig. [H]). However, these authors did not consider the radiative 
transfer in these structures and thus could not predict their optical properties 
like the Ha contrast against the disk. They only draw the black bars as other 
authors do in case of low-/3 dips. 



6.3. Other scenarios 

Finally, we should also mention that some authors considered waves propagating 
vertically within a flux-t ube and aimed at s upporting the cool plasma against 
the gravity. For example, iPecseli &: Engvoldl ((^00) have suggested MHD waves. 
They show an illustratory example (see their Fig. 1) of such "vertical flux-tubes" 
in a quiescent prominence with presumably vertical field in which the waves are 
pr opagating. H owever, this example rather resembles the situation discussed 
bv iLerovl (Il989l ). i.e., quasi- vertical fine-structure threads which are threaded 
by horizontal field lines. There is also a n ongoing debate whether the field 



i n barbs is more inclined or quasi- vertica l ( Zirker et al. 19981 ) or made of dips 



( Aulanier Sz Schmiedeil [20021 . Chae et al. 20051 ). In the former case, for fiows 



which do not correspond to a free fall, one would need a special support. 



7. Fine-Structure Energetics 

The energetics of quiescent prominences is still an unsolved problem. Although 
mutually connected, we can divide this problem in two aspects: the energetics 
of central cool cores of the FSE and the energetics of their PCTRs. 

The basic question is whether the incident radiation which almost fully de- 
termines the radiation properties of cool cores can also determine the kinetic 
temperature in these regions. In other words, is the cool core of FSE in a radia- 
tive equilibrium or do we need some extra heating to achieve the observed values 
of temperatures which are typically between 6000 to 9000 K. The computation 
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Figure 8. Gravity-induced magnetic dips as projecte d onto the disk. Eac h 
black bar corresponds to a vertical ID dip model. From lLow fc Petri^ (|2005f) . 



of the radiative equiHbrium under the non-LTE conditions is a rather difficult 
task, moreover, the result will strongly depend on the incident radiation fields 
used as boundary conditions for the solution of the radiation transfer problem. 
Radiation-equilibrium models in ID slabs, for a mixt i ire of hydrogen and helium 
plasma, were first constructed bv lHeaslv Sz Mihalas ( 19761 ) who arrived at rela- 
tively low temperatures, down to 4600 K. By adding some extra energy input via 
a heating function, higher temperatures were obtained in accordance with typical 
observations. Since the thermal conductivity is not efficient in these central cool 
parts, several other mechanisms were considered during the last few decades. 
Among them, we can mention mainly the wave dissipation or en thalpy trans- 
port. However, a surprisingly new result was recently obtained bv lGouttebrozd 
( 20071 ). who computed a set of fine-structure models in radiation equilibrium. 
ID axially-symmetric cylinders represent vertically standing thread s illuminated 
by th e disk radiation. Wit h currently used incident radiation fields ( Gouttebrozd 
20041 ) ■ ICouttebrozi ( 20071 ) arrived at r nuch higher radiati o n-equ ilibrium temper- 
atures as compared to the results of iHeaslv Sz Mihalad ( 19761 ) - we show this 
in Fig. O Namely thinner cylinders which are of interest for the fine-structure 
modeling (plots in Fig. [9^) reveal temperatures between 7000 and 9500 K for 
gas pressures between 0.5 to 0.01 dyn cm~^, respectively. This would then mean 
that no extra heating is required for cool central parts of fine-strucure threads. 
In these models only the radiative losses due to the hydrogen were considered. 



Fine Structure of Prominences 



283 




40 60 
Radius (km) 



1000 1500 
Radius (km) 



Figure 9. Temperature variation along the radius of vertical cylindrical 
threads in radiative equilibrium, for two diameters and various pressures. 
Pressures (in dyn cm~^) are indicated as curve labels. Left: for cylinders 
with t he width of 200 km; right: for a width of 5000 km. From [Gouttebroz^ 

(pool . 



which nevertheless can be compared with results of Heaslv &: MihalasI ( 19761 ) 
because the helium does not contribute much. However, there is still an open 
question concerning the importance of other spec ies like calcium, rnagnes ium or 
other optically-thin losses (see e.g., discussion by lAnze7 fc Heinzell ^199^ )). 

Concerning the energetics of the PCTR, our knowledge is also rather in- 
complete. At lower temperatures, say up to 10^ K, the c onduction along the 
magnetic field already plays a role ( Fontenla Sz Rovira 19841 ) . Moreover, in anal- 
ogy to CCTR ( Fontenla et al.|[l990l ). the ambipolar diffus i on wa s also considered 
for fine-structure prominence threads by iFontenla et all ( 19961 ). Quite recently, 
Vial et al.l (|2007l ) tested various kinds of ID models against the SOHO/SUMER 



observations of hydrogen Lyman line intensities and found that all models with 
the ambipolar diffusion predict unrealistically high integrated intensities of the 
Lyman /3 line. This can be related to th e relative orientatio n of the line of 
sight and the magnetic field as discussed in lHeinzel et al.l ( 20051 ). At still higher 
temperatures, most of the work was re stricted to the analysis of the differen- 



tial emission measure (DEM) - see e.g., Engvold ( 19981 ). However, this analysis 



is usually based on the so-called coronal approximation and thus the question 
arises whether the departures from this approximation, due to radiative exci- 
tations, can play any impo rtant ro l e. Fo r chromospheric conditions, this issue 
is discussed in this book by lAvrett ( 20071 ). At temperatures higher than 30000 
K, dynamical models of flux-tubes with prominence condensations were studied 
using extensive numer ical simulation s of th e time-dependent loop energetics and 
plasma dynamics (see Karpen et al. ( 20061 1 and references therein). 



8. RMHS Simulations 



The equations of magneto-hydrostatic equilibrium originally derived by KS be- 
come considerably simpler if one uses instead of the Cartesian coordinate x the 
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column-mass coordinate m, defined by the relation 

dm = —pdx , (2) 

where p is the the plasma density, m = at one surface of the prominence 
slab and m = M at the opposite. Loading of such mass into an initially hor- 
izontal magnetic field leads to a formation of a gravity-induced dip as we will 
demonstrate in the next Section. The pressure-balance equation for such a ID 
equilibrium dip configuration is governed by the equation 

in f Tfi \ 

pM=4pc^(l-^j+P0, (3) 

where p is the gas and turbulent pressure and pn is the coronal pressu re at 
the surfaces. This equation was first derived bv iHeaslv &: MihalasI ( 19761 ) and 



used to model prominences as a whole. At the slab surface one has the vertical 
component of the field vector = Bz\ which gives, together with the horizontal 
component B^ =const. 

M = . (4) 

2irg 

Using this formula, we obtain for pc 

Bl 2 " 8^ ■ 

The quantity pc can be interpreted in the following way: at the slab center we 
have the pressure 

Pccn =p(M/2) =Pc+PO- (6) 

If Po would be zero, then pcen = Pc = BIi/Stt, which is the magnetic pressure. 
Therefore, in this case the gas pressure at the slab centre will be equal to the 
magnetic pressure calculated with B = B^i- This formulation with the column 
mass has a great advantage of a simple analytical integration which is valid for 
any (i.e., non-constant) temperature and ionization-degree distribution. To get 
the density p{m) we use the state equation with the mean molecular mass 

l + 4a 

= T— — — rmn , (7 

1 + a + I 

where i is the ionization degree of hydrogen i = np/nn {np and nn are the 
proton and hydrogen densities, respectively), a the helium abundance relative 
to hydrogen and mn the hydrogen atom mass, i varies between zero (neutral 
gas) and unity (fully-ionized plasma). Inside the prominence and its PCTR, one 
can consider some schematic variation of i with depth but for a given prominence 
model the true ionization-degree structure results from rather complex non-LTE 

radiative-transfer calculations. 

The ID-slab MHS equilibrium of this type was first used bv lHeaslv Sz Mihalad 

( 19761 ). who combined it with the full set of non-LTE equations. Their models 



were aimed at describing the whole prominence. A similar study was repeated 
recently bv lAnzer Sz Heinzel ( 19991 ) who have considered a grid of models and 
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studied their energetics. In the latter work, new incident radiation fields and 
the partial-redist r ibutio n in hydrogen Lyman lines were used. As a next step, 
Heinzel Anzej ( 200ll ) have generalized these ID models to two dimensions 



and developed fully 2D MHS models coupled to the radiation field. The lat- 
te r were computed using the 2D radiative-transfer technique similar to that 
of [a uer k Paletoul (|l994l i In this way, all quantities including the ionization 
structure were consistently evaluated in the frame of such a 2D Radiation- MHS 
(RMHS) approach. However, contrary to previous work, these 2D models were 
aimed at representing the vertical fine-structure threads frequently observed in 
quiescent prominences. The only drawback of this kind of modeling is that 
Heinzel &: Anzer ( 200ll ) didn't consider the energy-balance problem and, in- 



stead, used some kind of empirical temperature structure. Their PCTR along 
the magnetic field lines is much more extended as compared to that across the 
lines, which corresponds to different thermal conductivities. Since these 2D 
models allow us to look at the fine-structure threads in various directions, we 
see that the respective synthetic spectra refiect very well the different structure 
of PCTR's. Namely the hydrogen Lyman lines appear quite reversed when we 
look across the field lines and much less reversed or even unreversed when look- 
in g along the field lines . This quite interesting behaviour was discovered already 



m 



Heinzel et all (|200lh on b asis of SOHQ/ SUMER spectra 



In subsequent papers, iHeinzel et al.l (|2005 ) and iGunar et al.l ( 20071 ) ex- 



tended this 2D modeling of vertical threads to a 12-level plus continuum hy- 
drogen model atom and studied in detail the formation of hydrogen Lyman lines 
and Lyman continuum. In Fig. [10] we show one of their results, and namely the 
2D contribution functions which illustrate the formation of hydrogen Lyman 
lines along and across the field lines, in a 2D vertical thread with PCTR tem- 
peratures ranging up to 50 000 K in this numerical box. This pattern is quite 
complex but leads to synthetic line profiles which are in good agreement with 
SOHO/SUMER observations (see also the paper by Gunar et al. on p. 317 ff 
in this volume). The sensitivity of Lyman-line profiles to orientation of the 
magnetic field w i th res pect to the line of sight has been also proven recently by 
Schmieder et al.l (|2007l ) who studied a kind of "round-shape filament" approach- 



ing the solar limb and consecutively showing its different parts above the limb 
having different orientations of the field lines. 

Finally, let us note that these 2D threads which can be anisotropically irra- 
diated represent the basic ingredient of a more realistic multi-thread modeling. 
Several such threads distributed in space will be illuminated by the solar ra- 
diation penetrating through this "forest" of threads and the mutual radiative 
interaction between all those threads can be consistently taken into account. 
Another important aspect of this 2D MHS modeling is that such threads or 
their clusters can represents prototypes of realistic models used to synthesize 
the Stokes profiles and thus study the infiuence of the fine structure on the 
polarization signals. 



9. RMHD Simulations 

Equation (3) gives the ID MHS equilibrium of the KS-type. In order to see how 
the cool and dense plasma structure can evolve in the magnetic field under the 
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Figure 10. Contribution functions along and across the field lines showing 
the formation depths of the hydrogen Lyman /3 (left two columns) and Lyman 
9 (right two columns). Two different models are displayed (twice three rows). 
In each part, the first row corresponds to the line center, the second one 
to the peak wavelength and the third one to the wing. Left columns refer 
to the direction along the field lines (x-direction) , while right columns to 
the direction across the field (^/-direction). The temperature contours are 
also drawn showing two quite different PCTR's. The line core of Lyman P 
is formed at the surface of the 2D thread (at the highest temperatures), the 
peaks deeper and in the line wing we can see the central parts of the structure. 
The Lyman 9 line is optically thinner and thus the thread more transparent in 
all wavelengths allowing the diagnostics of deeper layers. Note the different 
scales in both directions, in reality the "fibril" projected onto solar disk is 
much more stretched. 



action of Lorentz and gravity forces, Barta et aL ( 20071 ) solved numerically the 
system of compressible one-fluid MHD equations in a 2D vertical plane, start- 
ing with the dense and cool blob (a thread in 3D) representing a thread of the 
filament fine structure surrounded by the gravitationally stratified, constant- 
temperature hot corona. The initial magnetic field has only a horizontal com- 
ponent. At the very beginning the plasma blob starts to fall down, the internal 
electric current is induced inside it which generates the restoring Lorenz force. 
Due to the inertial mass the blob overshoots the (global) equilibrium point and 
an upward-directed Lorentz force prevails over the gravity turning the move- 
ment of the blob up. The system starts to oscillate but due to the numerical 
viscous term the oscillations are damped. After several oscillations the system 
is practically relaxed and in the MHS equilibrium. The evolution of the system 
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Figure 11. A sequence of four snapshots of the prominence fine-structure 
evolution according to the MHD system of equations. The simulation is car- 
ried out in the dimensionless units; lengths are expressed in the initial thread 
half- width Lt, times in the units of Alfven transit time ta = Lt/vA- Recalcu- 
lated for reasonable coronal conditions (coronal temperature Tc = 2 x lO^K, 
coronal gravitation length scale Lg = 6 x 10^ cm) we get Lt « 1200 km, 
TA ~ 20 s. Thus the thread characteristic dimension is ~ 2500 km and the 
final (relaxed) state in the last frame is reached after almost 20 hours. 



at four subsequent times is shown in Fig. [TTl These first simulations will now be 
replaced by fully 3D modeling which will also include the 3D radiative transfer. 
The latter is necessary to determine the ionization state of the plasma and its 
radiation losses. 



10. Comments on Magnetic Dips 

Any analysis of 2D images or movies can lead to misinterpretations just due to 
the lack of information on fully 3D structural and dynamical patter n. As an ex- 



ampl e we will mention here again the scenario of counter-streaming (jZirker et al 



19981 or Lin et al. 2003). A close inspection of high-resolution Ha movies (see 



them on CD enclosed to the Solar Physics Vol. 216) reveals that dark fine- 
strucure blobs (rather than long threads) move on relatively short paths, al- 
though sometimes in opposite directions. These motions don't resemble at all 
flows along long magnetic fluxtubes, but rather a kind of oscillations. Quite 
interestingly, when watching the TRACE movies of similar prominences, one 
clearly sees sideway oscillations of small blobs which, projected on the disk, 
would probably resemble the kind of motions seen on SST movies (projected 
blobs seen in Ha and TRACE 195 A are shown in Fig. [5]). However, these blobs 
visible on TRACE movies form very frequently the quasi-vertical threads also 
shown in Fig. [1] and Fig. [2] and thus, according to our discussion above, they 
should represent plasma condensations hanging in a dipped magnetic field which 
seems to oscillate in th e corona. This scenario is, however, quite different from 
that of iLin et al.l ( 20031 ) (flows along long fluxtubes without any dips). 



Another observational support for dips being located along the prominence 
spine is the fact that the Ha absorbing fibrils are concentrated along the spine 
and we typically don't see any significant absorption further away from the spine. 
This seems to indicate that there are no such flows which could be detected along 
the whole magnetic arcade which is of course much wider than the width of the 
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spine. Even Ha images or movies taken out of the line center don't show such 
a pattern. 

Actuahy there should be no principal controversy between low-/3 dip models 
and gravity-induced dip models. Yet there are no time-dependent simulations of 
how the plasma is condensed in a dipped field and how this may evolve leading to 
an enhancement of the plasma density and thus its weight which, eventually, will 
modify the magnetic structure of the dip. Existing time-dependent simulations 
of prominence condensations assume an initial shape of the magnetic loop and 
this is not changed during simulations even when the mass loading at the loop 
top increases. Static models of dips which neglect the prominence weight (fff 
extrapolations) can reasonably well reproduce the large-scale distribution of dips 
and this has proven to be a quite novel approach. On the other hand, static 
models of gravity-induced dips have been developed to understand the local 
magnetic topology of a dip and its thermodynamic and radiative properties. 
These two approaches should not be in conflict, but rather complementary: 
provided that the plasma-/? is large enough, the weight effects should be added 
to initially fff-dips. This is a challenge for future prominence modeling. As 
discussed above, today the actual values of plasma beta are still very uncertain 
and thus various configurations are possible. An important diagnostics test 
would be the quantitative modeling of the Ha contrast of fine-structure fibrils 
located along the filament spines. This contrast, which certainly increases with 
increasing spatial resolution of new telescopes, is due to the absorption of the 
background radiation in the Ha line (the scattering contribution to the source 
function is typically small). It has to be modeled in 3D fibril geometry which 
requires realistic estimate of t he mass loading. An exa mple of how this can be 
handled was given recently bv lHeinzel Sz Anzer ( 20061 ) . who used some kind of 
(2-|-l)D models to demonstrate the theoretical Ha contrast of gravity-induced 
magnetic dips as projected against the solar disk. A few examples are shown in 
Fig. [T2I where the fibril models do exhibit a significant stretching in the direction 
of the magnetic field and this resembles real dark fibrils as seen e.g., in Fig. [3l 



11. Conclusions 

To conclude this review, let us go again back to Secchi (1877). In his book Le 
Soleil one can read: "Les protuberances se presentent sous des aspects si bizarres 
et si capricieux qu'il est absolument impossible de les decrire avec quelque exac- 
titude. " After 130 years, this still seems to be the case. 
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